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ABSTRACT 

Recent  developments  in  pseudo-holography  provide  the  means  to 
generate  a  three-dimensional  image  of  an  object  from  any  radiation  it 
scatters.   One  promising  technique  uses  scattered  radiation  to  cast 
shadows  of  a  Fresnel  zone  plate,  referred  to  as  the  coded  image, 
which  are  recorded  on  film.   The  coded  image  is  then  decoded  by 
passing  coherent  light  through  it  which  is  diffracted  producing  an 
image  of  the  radiation  source. 

In  this  paper,  a  brief  illustration  of  pseudo-holography  using 
scattered  thermal  neutrons  is  presented,  along  with  a  development  of 
a  general  expression  for  the  coded  image.   The  remainder  of  the  paper 
is  devoted  to  developing  a  scattered  neutron  flux  source  function  by 
diffusion  theory  methods  for  a  non-reentrant,  steady  state,  homogeneous 
body  of  simple  geometry.   To  test  the  accuracy  of  the  source  function, 
the  computed  scattered  flux  is  compared  with  some  measured  values. 
The  comparison  supports  the  diffusion  approximation  for  detection 
methods  which  do  not  have  a  strong  forward  detecting  preference. 


I.   INTRODUCTION 

Recent  developments  in  pseudo-holography  provide  the  means  to 
generate  a  three-dimensional  image  of  an  object  from  any  radiation  it 
scatters.   One  promising  pseudo-holographic  technique  is  basically  a 
two-step  process.   The  scattered  radiation  is  used  to  cast  shadows  of 
a  Fresnel  zone  plate  (hereafter  termed  the  coded  image)  which  are 
recorded  on  film.   The  coded  image  is  then  decoded  by  passing  coherent 
light  through  it  which  is  diffracted  producing  an  image  of  the  radiation 
scattering  source.   This  technique,  first  proposed  for  radiation  sources 
by  Mertz  and  Young  in  x-ray  astronomy,  and  extended  to  gamma  ray 
sources  by  Barrett  '  '  '   and  Rogers  ,  has  been  successfully  applied  to 
actual  image  formation  using  x-rays  scattered  from  an  object  by  Sider, 

et  al.   The  use  of  scattered  neutrons  in  image  formation  is  presently 

g 
under  study  by  Rose,  et  al.   and  provides  the  motivation  for  this  paper. 

The  pseudo-holographic  technique  using  scattered  neutrons  differs 

from  the  use  of  scattered  x-rays  only  in  that  the  film  and  Fresnel  zone 

plate  respond  to  neutrons,  not  x-rays,  in  forming  the  coded  image.   It 

follows  that  an  analytical  representation  of  the  coded  image  can  be 

developed  by  describing  the  scattered  neutron  flux,  modifying  the 

scattered  flux  with  a  Fresnel  zone  plate  response,  and  applying  a  film 

(detector)  response  function  to  the  modified  flux  (shadows)  seen  by  the 

film.   The  first  step  in  this  analysis  is  an  accurate  description  of  the 

uninterrupted  scattered  flux.   It  is  the  intent  of  this  paper  to  develop 

an  analytical  method  which  describes  scattered  neutron  fluxes  emerging 

from  objects  in  a  beam  of  thermal  neutrons. 


First,  a  brief  illustration  of  pseudo-holography  is  presented 
along  with  a  development  of  a  general  expression  for  the  coded  image. 
The  remainder  of  the  paper  is  devoted  to  developing  a  scattered  flux 
source  function  by  diffusion  theory  methods  for  a  non-reentrant, 
steady  state,  homogeneous  body  of  simple  geometry.   To  test  the 
accuracy  of  the  source  function  the  computed  scattered  flux  is  compared 
with  some  measured  values. 


II.   BACKGROUND 

A.  Concept  Illustration 

The  basic  concept  of  the  pseudo-holographic  technique  noted  above, 
as  applied  to  thermal  neutrons,  is  illustrated  in  Figure  1.   A  Fresnel 
zone  plate  is  positioned  relative  to  the  scattering  body  and  a  suitable 
film-converter  camera  as  shown.   The  zone  plate  has  a  series  of  equi- 
area  annular  zones  which  are  alternately  transparent  and  opaque  to 
thermal  neutrons  and  consist,  for  example,  of  aluminum  and  cadmium, 
respectively.   Considering  only  a  point  scattering  source  of  thermal 
neutrons,  the  display  on  the  emulsion  film  (the  coded  image)  is  a 
shadow  of  the  zone  plate  the  shape  and  size  of  which  depend  monotoni- 
cally  on  the  position  of  the  source  in  the  z-plane,  (x,y);  the  position 
of  the  zone  plate  relative  to  the  source  plane,  z..-z;  and  the  position 
of  the  camera  relative  to  the  zone  plate,  z„-z  .   When  a  beam  of 
coherent  light  passes  through  the  coded  image,  the  light  is  diffracted 
into  a  focal  spot  which  is  an  image  of  the  original  point  scatterer. 
Considering  any  object  to  be  a  set  of  point  scatterers,  any  specimen 
can  be  imaged.   Each  plane  in  the  object  will  have  a  different  focal 
plane  when  imaged,  so  that  the  three-dimensional  information  in  the 
coded  image  can  be  viewed  tomographically. 

B.  Coded  Image 

Neglecting  the  detector  resolution*,  assuming  a  linear  film 
response,  and  noting  the  film  output  is  different  transparencies  for 


*The  detector  resolution  is  idealized  by  considering  the  film  and 
converter  to  be  in  the  same  plane  and  neglecting  the  film  resolution. 
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points  on  the  film,  the  coded  image  can  be  represented  by  the  neutron 

flux  shape  in  the  x"y"-plane.   As  shown  in  Figure  1,  (x",y")  describes 

points  in  the  film-converter  plane  at  z  =  z„;  (x',y')  describes  points 

in  the  zone  plate  plane  at  z  =  z  ;  and  (x,y,z)  describes  points  in 

the  scattering  body. 

Let  F  (x,y,z,£2),  the  scattering  rate  density,  represent  the 
s       — 

number  of  neutrons  per  unit  time  per  unit  volume  per  solid  angle 
scattered  in  the  volume  dxdydz  about  x,y,z  and  in  the  direction  dJ] 
about  fl.   Let  G(x,y  ,z,£2  -*■  x"  ,y"  ,z_  ,£3")  represent  the  angular  flux 
Green's  function  to  the  film-converter  plane,  then 


-ri(r,r") 

~  -         r"-r 


G(x,y,z,ft  +  xM,y",z9,n")  =  ^ 5-  6 

1  |r-r"| 


n  - 1=~7 

—    r-r 


6  (£"-£)    (i) 


where  _r,r_'  and  r_"  are  vectors  from  the  origin  (see  Figure  1)  to  points 
in  the  scattering  body  (x,y,z),  to  points  in  the  zone  plate  plane 
(x',y',z'),  and  to  points  in  the  film-converter  plane  (x",y",z„), 
respectively;  and  r)(r_,r_")  ^s  tne  "optical  path  length"  between  r_  and  r_". 

In  this  case  it  is  convenient  to  consider  l")(r,r")  in  two  parts, 
the  "optical  path  length"  in  the  scattering  body,  ru  (r_,£M)  ,  and  the 
"optical  path  length"  outside  the  body,  r\    (r_,_r").  Then 


-n(r,r")    -n  (r,r")    -n9(r,r") 
;   —  —   —  e  x  —  —   •  e  i.   —  — 


(2) 


where  the  second  exponential  on  the  right  hand  side  represents  the  zone 
plate  transparency  which  is  defined  by  a  set  of  step  functions  or  a 
series  expansion  . 


The  neutron  flux,  $(x" ,y" ,z„) ,  at  the  film-converter  plane  is 

then 

(3) 


$(X",y,,,z2)  = 


F  (x,y,z,ft)G(x,y,z,ft  -*■   x"  ,y",z0£")dfidg"dr 
S  -  2. 


Q   ft"Vol 


If  an  angular  emitted  neutron  current  from  the  surface  of  the  body  can 

be  described  as  K  (x,y,z,fl),  then 

S  (4) 


$(x",y",z2)= 


K  (x,y,z,n)G(x,y,z,£  -»■  x",y",z2,fi")d>gclf2!ldr 


where 


Q  Q"   S 


represents  integration  over  the  body  surface,  dr  is  a  surface 


s 


element  of  the  body  and  in  this  case  n,  (r,r_")  =  0.   Note  that  with 
the  zone  plate  removed,  ri  (r_,r_")  =  0  in  equations  (3)  and  (4)  above, 
and  the  uninterrupted  scattered  neutron  flux  at  (x",y")  in  the  film- 
converter  plane  is  described. 

To  evaluate  equations  (3)  or  (4)  an  expression  for  the  emitted 
neutron  current  or  the  scattered  neutron  source  must  be  developed  via 
a  suitable  neutron  transport  approximation.   The  accuracy  of  the  source 
function  can  be  verified  by  computing  the  uninterrupted  scattered  flux 
from  equations  (3)  or  (4)  and  comparing  its  shape  with  measured  values. 
As  noted  above,  the  flux  shape  in  the  (x",y")  plane  is  sufficient  to 
represent  the  coded  image. 


III.   SCATTERED  FLUX  FROM  A  REGULAR  SOLID 

For  the  development  of  a  coded  image  of  the  interior  structure  to 
be  practicable,  the  thermal  neutron  scattering  cross  section,  £  ,  of  the 
scattering  object  must  be  much  greater  than  the  thermal  neutron  absorp- 
tion cross  section,  £  .   This  suggests  the  possible  application  of  the 
diffusion  approximation  to  neutron  transport  in  developing  a  scattered 
source  function  as  used  in  equations  (3)  or  (4) ,  where  some  approximation 
will  be  required  to  define  the  directional  dependence  of  the  scattered 
source  function. 

A.   Scattered  Source  Function 

It  is  assumed  that  an  adequate  test  case  is  to  examine  the 

scattered  flux  from  a  regular  solid  in  a  beam  of  thermal  neutrons 

of  intensity  I   (neutrons  per  unit  time  per  unit  area)  as  shown  in 

Figure  2,  where  the  composition  of  the  solid  is  homogeneous  and  such 

that  L  >>  £  .   As  mentioned  earlier,  either  a  scattering  source  function 
s    a 

or  an  angular  emitted  current  source  function  can  be  used  in  computing 

the  scattered  flux.   To  develop  a  scattering  source  function,  the  flux 

in  the  solid,  $(x,y  tz) ,      is  derived  by  diffusion  theory  methods,  and 

the  solid  is  considered  an  isotropic  volumetric  source  of  strength 
ig$(x,y,z) 


4tt      .   The  scattered  flux  at  some  point  away  from  the  solid, 

(x  ,y  ,z  )  in  Figure  2  would  be  the  uncollided  contribution  from  each 
ooo 

point  in  the  solid.   An  angular  emitted  current  source  function  on  the 
surface  of  the  solid  is  defined  directly  from  the  emitted  current 
derived  via  diffusion  theory,  where  some  emission  angular  dependence 
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must  be  assumed.   The  use  of  diffusion  theory  in  computing  reflection 
coefficients  and  the  relative  ease  in  computing  the  scattered  flux  from 
an  angular  emitted  current  source  suggests  the  use  of  the  angular 
emitted  source. 

The  diffusion  equation  to  be  solved  is 


DV  $(x,y,z)  -  £  *(x,y,z)  +  S(x,y,z)  =  0 


(5) 


where  D  is  the  diffusion  coefficient  of  the  solid  and  S(x,y,z)  is  an 
isotropic  volumetric  source.   It  is  convenient,  in  terms  of  describing 
suitable  boundary  conditions,  to  let  S(x,y,z)  represent  the  first 
collision  scattering  density,  i.e., 


-E  z 
S(x,y,z)  =  I  I   e   t 
s  o 


(6) 


where  £  is  the  total  thermal  neutron  cross  section  of  the  solid  and  it 
is  assumed  that  the  first  collision  scattering  density  is  isotropic. 
The  boundary  conditions  are  that  $(x,y,z)  is  zero  at  the  extrapolated 
boundaries  of  the  solid. 

The  solution  to  equation  (5)  (see  Appendix  A)  obtained  by  the 
method  of  eigenfunctions  for  the  geometry  in  Figure  2  is 


>(x,y,z)  =  I     A 

lmn 


lmn 


sin 


Ittx 
c 

sin 

mTTy 

sin 

nirz 

a 

fe 

(7) 


10 


where 


32   Z   I     K 
s   o 


lmn 


i.a«A£_,   B-^m]2 


K  -   *-hb 


=        .      rim(B-2D)l        nil  [i 

-  lt  sin  [ b — J  "  T   cos  |_~ 


,    s         .      rnTT2D1     ,    niT  pttr2Dn 

+  lt  sin  l__j  +  _  cos  |_-^g — J 


L  =  D/Z    ,    B   =  b+4D,      c    =  c+4D,    I  =   a+4D 
a 


T_mn 


■m 


mTT 


nil 


t   2 


Noting  the  symmetry  in  Figure  2,  x  =  c/2  and  y  =  a/2  is  the  center 

of  symmetry  in  the  xy-plane,  and  considering  points  (x  ,y  ,  z  )  which  do 

not  see  the  ends  of  the  solid,  i.e.,  2D<z  <b+2D,  then  the  emitted 

o 

current  need  only  be  defined  for  faces  1  and  2  of  the  solid.  The  emitted 
current  from  face  1,  J   (c+2D,y,z),  by  Fick's  Law  is 


J+  (c+2D,y,z)  -  »(c+2D,y,z)  -f 


9$(x,y,z) 

8x 


(8) 


x  =  c+2D 


where 


3$(x,y,z) 
3x 


x  =  c+2D   lmn 


lT7_ 
C 


A,    cos 
lmn 


1tt(c+2D) 
c 


sin 


mTry 
a 


sin 


nTTZ~l 

_b  J 


11 


and  the  emitted  current  from  face  2  is 


J   (x,a+2D,z) 


>(x,a+2D,z)    D 
4       "2 


5$(x,y,z) 
3y 


J  y  =  a+2D 


(9) 


where 


9$(x,y,z) 
3y 


I     -w—  A,    cos 


J  y  =  a+2D 


lmn 


a   lmn 


mTT(a+2D) 
a 


]  -  m 


rmTzi 
sin  |— g— J 


B.   The  Scattered  Flux 

To  use  equation  (4),  with  ri(r_,r_")  =  0,  in  calculating  the 

uninterrupted  scattered  flux,  $(x  ,y  ,z  )  at  (x  ,y  ,z  ),  it  is  necessary 

o  o  o       o  o  o 

to  define  some  angular  dependence  to  the  emitted  current.   Realizing 
it  is  difficult  to  make  any  generalization  about  the  angular  dependence 
of  the  emitted  current,  two  possibilities  for  defining  an  angular 
dependence  in  terms  of  isotropy  are: 

1.  Replace  the  emitting  face  with  a  plane  semi-isotropic  source 
(hereafter  termed  the  semi-isotropic  emitted  current  assumption)  of 
strength  J  (x,y,z)  defined  by  equation  (8)  or  (9). 

2.  Assume  the  current  density,  j(x,y,z,£2),  is  semi-isotropic 
(hereafter  termed  the  semi-isotropic  current  density  assumption) 
then 


J   (x,y,z)  = 


n.   j(x,y,z,fi)dfi 


(10) 


2tt 


12 

where  j(x,y,z,£2)  =  Cfi  and  11  is  the  unit  vector  normal  to  the  emitting 
area. 

Based  on  the  semi-isotropic  emitted  current  assumption  the  angular 
emitted  current  from  face  1,  K  (c+2D,y,z,£jO  ,  is 

S  J- 

J+(c+2D,y,z) 
K   (c+2D,y,z,fi)  - 


slv   "'  '-        2tt  (11) 


and  the  angular  emitted  current  from  face  2,  K  (x,a+2D,zfi),  is 


J+(x,a+2D,z) 
K  .(x,a+2D,z,n)  =  -2 — (12) 

Si.  —  Alt 


Based  on  the  semi-isotropic  current  density  assumption,  the 
angular  emitted  current  from  face  1,  K  „ (c+2D,y ,z,fi!)}  is 

K  ,  (c+2D,y,z,ft)  =  n.  j  (c+2D,y  ,z  ,Q)    =  C.  cos  d>  (13) 

si         —    —  —         —     1 

where  cos  <j>  ■  il.  fi  and  §   is  then  the  angle  between  the  neutron  direction 
and  unit  normal  to  the  emitting  face.   Evaluation  of  equation  (10) 
yields 

Ac+2D,y,z) 
C,  = 


1        it  (14) 


Similarly  the  angular  emitted  current  from  face  2,  K   (x,a+2D,z ,fi) ,  is 

s  ^ 

J+(x,a+2D  z) 

K  .  (x,a+2D,z  ,Ci)    =  — * cos  d)  ,--» 

s2         —         u  (15) 


13 


Using  equation  (4),  with  r)(r_,r")  =  0,  the  uninterrupted  flux, 

$(x  ,y  ,z  )  at  (x  ,y  ,z  ),  recalling  that  2D<z  <b+2D,  is  determined 
ooo       ooo  o 

for  the  two  angular  emitted  currents  as  follows : 

1.  For  the  semi-isotropic  emitted  current  assumption: 

a.   If  2D<y  <a+2D  and  x  >c+2D  then  only  face  1  contributes 

J  o  o  ' 

to  the  flux  at  (x  ,y  ,z  )  and  evaluation  of  equation  (4)  over 
ooo 


Q  and  £2"  yields 


^wV  = 


a+2D  b+2D  , 

J    (c+2D,y,z)dydz 


2D   2D      2TT  r[xo-(c+2D)]2+(yo-y)2+(zo-z)2~| 


b.   Similarly,  if  2D<x  <c+2D  and  y  >a-t-2D  then  only  face  2 

o  o 

contributes  to  the  flux  at  (x  ,y  ,z  )  and 

o  o  o 


$(xo'yo'Zo}  = 


C+2D  b+2D        + 

J  (x,a+2D,z)dydz 


(16) 


(17) 


2D  2D  2Tr|"(x  -x)2+[y  -(a+2D)]  2+(z  -z)2] 


c.   If  y  >a+2D  and  x  >c+2D,  then  faces  1  and  2  contribute 
o  o 

to  the  flux  at  (x  ,y  ,z  )  and 
ooo 


"Kx  ,y  ,z  )  =  equation  (17)  +  equation  (16)  (18) 

ooo 


2.   For  the  semi-isotropic  current  density  assumption: 

a.   If  only  face  1  contributes  to  the  flux  at  (x  ,y  ,z  )  the 

ooo 

evaluation  of   equations    (4)    over  Si  and  fi"   yields 


14 


a+2D  b+2D 


$(x  ,y  ,z  )  = 
o  o  o 


cos  (J)  J  (c+2D,y,z)dydz 
2D  2D   7Tr[xo-(c+2D)]  2+(yo-y)2+(zo-z)2l 


(19) 


where  the  integration  of  equation  (4)  over  &   requires  that 


cos  $..  =  cos  cj> 


r"_r 
-"  |r"-r| 


x  -(c+2D) 
o 


r[xo-(c+2D)]  2+(yo-y)2+(zo-z)2l 


2-|  1/2 


b.   If  only  face  2  contributes  to  the  flux  at  (x  ,y  ,z  )  then 

o  o  o 


c+2D  b+2D 


,(xo'yo'Zo)  = 


cos  (j)„J  (x,a+2D,z)dydz 
2D   2D  tt  f(x  -x)2+[yQ-(a+2D)]  2+(z  -z)2"| 


(20) 


where 


cos  <j>„  = 


y  -(a+2D) 
o 


Rx  -x)  +fy  -(a+2D)l  +(z  -z)2 
I   o     L  o      J     o 


1/2 


c.   If  faces  1  and  2  contribute  to  the  flux  at  (x  ,y  ,z  )  then 

o  o  o 


>(x  ,y  ,z  )=  equation  (19)  +  equation  (20) 
o  o  o 


(21) 


To  evaluate  the  above  integrals,  a  numerical  technique  is 
conveniently  applied.   The  series  in  equations  (8)  and  (9)  are  uniformly 
convergent  so  that  the  series  can  be  integrated  term  by  term.   Thus  each 
term  in  equations  (16)  ,  (17) ,  (19)  and  (20)  can  be  integrated  numerically 


15 

and  a  suitable  convergence  established.   Results  of  scattered  fluxes 
for  a  specific  regular  solid  computed  as  outlined  above  are  contained 
in  the  next  section. 
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IV.   METHOD  OF  MEASUREMENT  AND  CALCULATION 

A.   Method  of  Measurement 

Measurements  of  the  scattered  flux  were  obtained  with  an  arrangement 

simulating  the  geometry  of  Figure  2.   As  shown  in  Figure  3,  a  regular 

solid  specimen,  a  polyethylene  (C  H.)  block  (2.86  cm  square  x  5.16  cm 

long),  was  placed  in  a  beam  of  thermal  neutrons.   A  BF  detector  was 

used  to  count  scattered  neutrons.   The  detector  longitudinal  axis 

remained  perpendicular  to  face  2  and  in  the  yz-plane  at  x  =  1.43  cm  + 

2D,  where  D  is  the  diffusion  coefficient  in  cm.   The  center  of  the 

detector  window  was  located  at  the  following  points  (x  ,y  ,z  )  for 

o  o  o 

measurement : 

1.  (1.43  cm  +  2D,   7.86  cm  +  2D,z  )  where  z  was  incremented  by 

o        o 

by  .5  cm  from  .5  cm  +  2D  to  4.5  cm  +  2D. 

2.  Same  as  1,  except  y  =  12.86  cm  +  2D. 

The  Pennsylvania  State  University  TRIGA  Reactor  (Breazeale 
Nuclear  Reactor  Facility)  provided  the  source  neutrons.   A  D„0  tank 
between  the  beam  port  entrance  and  the  reactor  core  provided  sufficient 

thermalization  such  that  the  thermal  flux  at  the  beam  port  exit  at  a 

7  2 

reactor  power  of  100KW  was  on  the  order  of  10  neutrons  per  cm  per 

second.   Collimation  was  achieved  by  placing  in  the  beam  port  two  iron 

plugs  with  1-1/2  inch  holes  followed  by  a  copper  plug  with  a  3/4  inch 

hole.   Smaller  collimations  would  reduce  background  scattering,  however, 

to  maintain  the  specimen  totally  within  the  beam  the  above  sizes  were 

the  minimum  allowable.   The  specimen  was  placed  15  inches  from  the  beam 

port  exit  and  at  this  point,  it  was  estimated  that  the  beam  had  diverged 
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to  approximately  twice  the  exit  diameter  with  a  flux  on  the  order  of 
10  neutrons  per  cm  per  second.   The  cadmium  ratio  in  the  beam  was 
approximately  50. 

The  BF  detector  was  14.5  inches  long  with  a  1  inch  outer  diameter 
and  was  filled  to  a  pressure  of  60  cm  Hg  with  BF  96%  enriched  in  B 
The  detector  was  shielded  with  a  cadmium  jacket,  and  the  window  had  a 
cadmium  cover  with  a  1/16  inch  diameter  hole.   By  placing  a  neutron 
shutter  between  the  beam  port  exit  and  the  detector,  an  average  signal 
to  noise  ratio  of  2  was  achieved.   As  shown  in  Figure  3  the  neutron 
counting  arrangement  was  typical  for  a  gas  filled  proportional  counter. 

The  polyethylene  specimen  was  held  in  place  by  an  aluminum  clamp 
attached  to  a  horizontal  bar  1/4  inches  in  diameter.   Three  of  the 
side  faces  were  covered  with  a  cadmium  sheet  which  reduced  background 
scatter  from  the  aluminum  holding  stand,  and  maintained  the  integrity 
of  the  non-reentrant  boundaries. 

One  minute  background  counts  were  taken  at  each  point  described 

above  with  only  the  specimen  removed.   The  gross  counts  were  also 

taken  for  one  minute  and  were  sufficiently  high  to  yield  counting 

fractional  standard  deviations  on  the  order  of  .025.   The  counts 

obtained  for  points  5  cm  (y  =  7.86  cm  +  2D)  and  10  cm  (y  +  12.86  cm  + 

2D)  from  face  2  were  normalized  with  the  z   -  2D  =  2.5  cm  point.   The 

o 

results  are  displayed  on  Figures  4  and  5.   As  discussed  earlier,  the 
normalized  curve  represents  the  relative  scattered  flux  shape  as 
amended  by  the  detector  response. 

Additional  measurements  were  made  at  2.5  cm  (y  =  5.36  cm  +  2D) 
from  face  2.  The  measurements  were  repeated  with  various  strips  of 
cadmium  on  face  2  as  follows: 
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1.  A  1  cm  wide  strip  beginning  at  x  =  1.5  cm  +  2D  and  ending  at 
x  =  2.5  cm  +  2D. 

2.  A  2  cm  wide  strip  beginning  at  x  =  1.5  cm  +  2D  and  ending  at 
x  ■  3.5  cm  +  2D. 

All  counts  obtained  were  normalized  with  count  obtained  at  (1.43  cm 
+  2D,  5.36  cm  +  2D,  2.5  cm  +  2D).   The  results  are  displayed  on  Figures 
6 ,  7 ,  and  8 . 

B.   Method  of  Calculation 

The  method  of  calculation  was  outlined  in  the  previous  section, 
and  in  this  case  equations  (17)  and  (20)  applied.   To  compare  the 
calculations  with  the  measured  values,  a  response  function  was  inserted 
representing  the  BF  detector  response.   The  response  function  was 
determined  as  follows. 

Neglecting  the  finite  area  of  the  detector  face,  the  probability 
of  detection  for  a  neutron  entering  the  detector  at  an  angle  0  with 
the  detector  face  is 

1  -  exp(-Z  r/cos  6)  (22) 

-B   -B  10 

where  Z  =  £  is  the  thermal  neutron  absorption  cross  section  of  B 
t    a  r 

at  60  cm  Hg,  r  is  the  radius  of  the  detector  and  r/cos  6  represents 
the  path  length  of  the  neutron  in  the  detector.   For  the  geometry  in 
Figure  3 


cos  0  = 


2       2        -,1/2 

(x  -x)  +  (z  -z) 
o       o 


(x  -x)  +  [y  -(a+2D)l  +  (z  -z) 
o      L:  o      Jo     J 


(23) 
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Considering  the  finite  area  of  the  detector  window,  a  unit  flux 
in  the  6  direction  to  the  detector  face  produces  a  current  of  sin  9 
entering  the  detector,  where  from  equation  (23) 


sin 


[y0-(a+2D>] 


(x  -x)2+ry  -(a+2D)l2+(z   -z)2 
o  L  o  J  o 


1/2 


(24) 


The  detector  response  function,  R(x  ,y  ,z  ,x,y,z),  is  then  the  fraction 
of  the  flux  seen,  equation  (24),  times  the  probability  of  detection, 
equation  (22) . 

Assuming  a  semi-isotropic  emitted  current  and  accounting  for  the 
detector  response,  equation  (17)  was  modified  to 


c+2D  b+2D 


$(x  ,y  ,z  )  = 
o  o  o 


R(x  ,y  ,z  ,x,y,z)  J  (x,a+2D, z)dxdz 
o   o   o   J  y   ' 

2D  2D  27Tr(xo-x)2+[^o-(a+2Dyf+(zo-z)2  "| 


(25) 


Assuming  a  semi-isotropic  current  density  equation  (20)  was  modified  to 


c+2D  b+2D 


(26) 


>(x    ,y    ,z   ) 
o     o      o 


R(x    ,y    ,z    ,x,y,z)cos    6„   J    (x,a+2D,z)dxdz 
ooo z     y 

2D   2D      tt  |(x  -x)2+[y   -(a+2D)]2+(z   -z)2] 


Equations  (25)  and  (26)  were  solved,  as  noted  previously,  by  integrating 
each  term  numerically  using  as  IBM  system  370  Mod  168  computer  and 
Argonne  National  Laboratory  subroutine  DROMB   .   The  series  was  truncated 
when  the  higher  mode  integral  terms  were  less  than  .1%  of  the  existing 
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sum.   The  thermal  neutron  diffusion  properties  of  polyethylene  were 
obtained  from  the  ANL-5800   .   Calculations  were  made  for  each  point 
measured  by  the  detector  and  were  similarly  normalized.   The  results 
for  the  semi-isotropic  emitted  current  and  for  the  semi-isotropic 
current  density  are  displayed  on  Figures  4  and  5. 

The  calculations  were  repeated  for  the  detector  at  2.5  cm  from 
face  2  as  used  in  obtaining  the  measured  values.   Calculations  were 
made  at  each  point  measured  by  the  detector  and  were  similarly 
normalized.   Calculations  with  the  cadmium  strips  on  face  2  were  made 
the  semi-isotropic  current  density  assumption  only.   The  results  are 
displayed  on  Figures  6,  7,  and  8. 
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V.   ANALYSIS  OF  RESULTS 

In  analyzing  the  results  on  Figures  4  through  8  it  is  helpful  to 
bear  in  mind  that: 

1.  The  detector  had  a  forward  counting  preference  which  decreased 
as  the  detector  was  moved  away  from  the  emitting  surface  (see  equations 
(22)  and  (24)). 

2.  The  emitted  current,  used  in  determining  the  angular  emitted 
current,  resulted  from  a  diffusion  calculation. 

3.  The  angular  emitted  current  was  determined  from  assumptions 
regarding  angular  dependence. 

As  shown  by  Figure  4,  at  10  cm  from  the  emitting  face  the  measured 
flux  shape  was  generally  well  represented  by  the  calculated  flux  shapes 
with  differences  only  at  points  near  the  boundaries  of  the  emitting 
face.   Coupled  with  the  detector  forward  counting  preference,  the 
boundary  errors  would  be  expected  from  a  diffusion  calculation.   The 
differences,  as  displayed  in  Figure  4,  indicate  that  the  boundary 
errors  were  the  result  of  calculated  emitted  currents  which  were  low 
near  the  edges  of  the  emitting  face. 

Additionally,  a  portion  of  the  observed  differences  can  be 

attributed  to  the  assumptions  used  in  defining  an  angular  emitted 

current  as  the  calculated  flux  shape  for  the  semi-isotropic  current 

emission  differs  from  the  calculated  flux  shape  for  the  semi-isotropic 

current  density  emission.   To  quantitate  somewhat,  the  maximum 

difference  between  the  calculated  relative  flux  for  a  semi-isotropic 

emitted  current  and  the  measured  relative  flux  occurred  at  z  -  2D  = 

o 

4.5  cm  and  was  15.1%.   The  difference  between  the  calculated  relative 
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flux  for  a  semi-isotropic  current  density  and  the  measured  relative 
flux  at  the  same  point  was  18.3%.   Thus,  the  errors  displayed  by 
Figure  4  were  due  primarily  to  errors  in  the  calculated  emitted  current 
near  the  boundaries . 

As  shown  by  Figure  5,  when  the  detector  was  5  cm  from  the  emitting 
face,  the  boundary  errors  increased.   This  would  be  expected  as  the 
detector  preference  for  points  directly  in  front  increased  and  the 
calculated  emitted  current  was  in  error  at  the  boundaries.   Moreover, 
errors  in  the  assumed  angular  dependence  of  the  emitted  current 
increased.   To  quantitate  again  at  z  -  2D  =  4.5  cm,  the  difference 
between  the  calculated  relative  flux  for  a  semi-isotropic  emitted 
current  and  the  measured  relative  flux  was  26.8%.   The  difference 
between  the  calculated  relative  flux  for  a  semi-isotropic  current 
density  and  the  measured  relative  flux  was  33.8%. 

As  shown  by  Figure  6,  when  the  detector  was  2.5  cm  from  the 
emitting  face,  the  differences  increased  further  as  would  be  expected. 
The  differences  in  the  calculated  values  for  the  semi-isotropic  current 
density  and  the  semi-isotropic  current  were  small  compared  to  the  total 
errors.   The  principal  source  of  the  errors  was  again  the  calculated 
emitted  current.   Moreover,  it  should  be  noted  that  the  block  dimensions 
were  only  1L  x  2L,  where  L  represents  the  diffusion  length. 

The  measured  and  calculated  flux  shapes  at  2.5  cm  from  the  emitting 
face  with  1  cm  and  2  cm  cadmium  strips  on  the  emitting  face  are 
displayed  on  Figures  7  and  8,  respectively.   The  boundary  errors  dis- 
played by  Figure  6  are  not  as  evident  here,  due  to  the  normalization. 
In  this  case  a  significant  portion  of  the  face  which  had  a  reasonably 
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correct  calculated  emitted  current  was  eliminated  by  the  cadmium  strip. 

As  a  result,  the  calculated  flux  at  the  z  -  2D  =  2.5  cm  point  had  a 

o  r 

larger  contribution  from  the  boundary  areas  yielding  a  low  value. 
When  normalized  to  this  low  value,  the  low  values  near  the  edges  of 
the  emitting  face  were  dampened,  decreasing  the  errors  in  the 
calculated  relative  flux. 

The  above  comments  can  be  summarized  as  follows.   Differences 
between  the  calculated  values  and  the  measured  values  were  due  primarily 
to  boundary  errors  in  the  calculated  emitted  current,  with  very 
minor  effects  due  to  the  arbitrary  assumptions  used  in  defining  the 
angular  dependence  to  the  emitted  current.   As  the  detector  forward 
preference  was  decreased,  the  observed  differences  decreased  and  the 
calculated  flux  shapes  represented  well  the  measured  values.   Although 
large  differences  existed  between  the  calculated  and  measured  values 
when  the  detector  had  a  strong  forward  preference,  measured  dips  in  the 
relative  flux  over  a  strong  absorbing  strip  were  represented  by  the 
calculated  values. 
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VI .      CONCLUSIONS 

A.  Summary 

The  diffusion  approximation  was  applied  in  calculating  emitted 
neutron  currents  from  a  homogeneous  regular  solid  in  a  beam  of  thermal 
neutrons.   The  calculated  emitted  current,  with  an  assumed  angular 
dependence,  was  used  in  computing  the  scattered  flux  shape  about  the 
solid.   Calculated  results,  modified  for  the  response  of  a  cylindrical 
gas  filled  proportional  counter,  were  compared  to  values  measured  with 
a  BF  proportional  counter  at  points  which  received  flux  contributions 
from  only  one  face  of  the  solid.   The  observed  differences  in  the 
calculated  and  measured  flux  shapes  were  due  primarily  to  errors  in 
the  calculated  emitted  current  about  the  edges  of  the  face,  with  very 
minor  effects  due  to  the  assumed  angular  dependence  of  the  emitted 
current.   As  the  directional  preference  of  the  detector  decreased, 
the  observed  differences  decreased.   Moreover,  regardless  of  a  signifi- 
cant directional  detector  preference,  trends  in  flux  shapes  over  strong 
absorbing  strips  were  well  represented.   It  can  be  concluded  that  for  a 
detector  with  a  lesser  directional  preference,  such  as  a  film-converter 
detector,  the  calculated  flux  shape  would  better  represent  the 
measured  flux  and  would  be  sufficient  to  complete  the  analytical 
representation. 

B.  Suggestions  for  Future  Work 

The  description  of  the  scattered  neutron  flux  shape  represents 
the  first  step  in  an  analytical  representation  of  the  pseudo-hologram. 
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Future  work  is  suggested  to  complete  the  analytical  representation  by 
including  the  response  functions  for  the  Fresnel  zone  plate  previously 
described  and  for  a  typical  film-converter  detector.   The  detector 
response  function  will  have  to  be  suitably  approximated  as  was  done 
herein  for  the  BF  detector  response  function.   The  flux  shapes  so 
determined  could  then  be  compared  to  calculated  flux  shapes  based  on 
an  ideal  detector  and  an  ideal  Fresnel  zone  plate  in  an  effort  to 
develop  a  correction  factor  for  pseudo-holograms  practically  obtained. 

One  possible  application  of  a  correction  factor  would  be  in  correcting 

12 

pseudo-holograms  which  are  decoded  numerically  on  a  computer  system 

The  correction  factor,  depending  on  a  general  application,  would  at 
least  decrease  some  of  the  resolution  loss  experienced  in  the  real 
system. 
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APPENDIX 
SOLUTION  OF  DIFFUSION  EQUATION 

The  geometry  is  as  shown  in  Figure  2.   Then  the  extrapolated 
boundaries  are 

c  =  c  +  4D 
a  =  a  +  4D 
B  =  b  +  4D 
where  c,  a  and  b  are  the  block  dimensions. 
The  equation  to  be  solved  is 


DV2$(x,y,z)  -  S  *(x,y,z)  +  S(x,y,z)  =  0  (A-l) 

3. 


where  the  source  is 


S(x,y,z)  =   / 


so  for  2D<z<b-2D  (A-2) 


0     otherwise 


and  the  boundary  conditions  are  $(x,y,z)  extrapolated   -  ,.    _x 

boundary 

To  solve  (A-l)  use 


V2c)>(x,y,z)  +  B^cKx.y.z)  =  0  (A-4) 


to  generate  a  linear  independent,  complete,  orthogonal  set  of  functions 
{<f>.   (x,y,z)}  for  l,m,n  =  l,2,3,-««-  where 
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.     ,1ttx.     .     /miry\    .     /nTTzN 
^Imn  =   ■la(-T)«l»(-r)8ln(T") 


and 


(A-5) 


lmn 


(ll)2  +   (™L)2  +   (-1)2 
cab 


(A-6) 


Then 


S(x,y,z)   =   I      S.        4>        (x.y.z)        and 


lmn 


lmn     lmn 


(A-7) 


*(x,y,z)    =   I     A.         $         (x.y.z) 


lmn 


lmn     lmn 


(A-8) 


From  the  orthogonality  property  of  {<)).   (x,y,z)}  and  equations 
(A-2)  and (A-7),  then 
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s  o 


lmn   s-,   2  rv2    .nTT  2l 
blmTT  [Zt  +  (g-)  j 


(A-9) 


where 


K  =  e 


-I,b 


■Z  sin  [ 


TTn(E-2D) 


nTT 


:os 


mr(E-2D) 


+ 


a   .  ,nTT2D,  ,  nTT     /niT2D- 
Ztsxn(^— )  +rcos(-r-) 


(A-10) 


Substituting  equations  (A-7)  and  (A-8)  into  equation  (A-l) ,  using 
equation  (A-4)  and  the  linear  independence  of  {((>..   (x,y,z)},  then 


lmn 


lmn 


2  2 
Z    (1+L  B7   ) 
a      lmn 


(A-ll) 
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where 


L2  =  D/Z 
a 


Then  using  equations  (A-8) ,  (A-9) ,  (A-10)  and  (A-ll) 


$(x'y'z)  *l     Almn  W*'y'z)  (A"12) 

lmn 


35 


BIBLIOGRAPHY 


1.  Mertz,  L.  and  N.  0.  Young.   "Fresnel  Transformation  of  Images." 
Proceedings  of  the  International  Conference  on  Optical 
Instrumentation,  London,  305  (1961). 

2.  Barrett,  H.  H.   "Fresnel  Zone  Plate  Imaging  in  Nuclear  Medicine." 
Journal  of  Nuclear  Medicine,  13:382  (1972). 

3.  Barrett,  H.  H. ,  P.  T.  Wilson  and  G.  D.  DeMeester.   "The  Use  of 
Halftone  Screens  in  Fresnel  Zone  Plate  Imaging  of  Incoherent 
Sources."  Optics  Communications  5:398  (1972). 

4.  Barrett,  H.  H.  and  F.  A.  Harrigan.   "Theory  of  Fresnel  Zone  Plate 
Imaging  of  Gamma  Rays."   Raytheon  Technical  Report  T926. 
Submitted  to  Applied  Optics. 

5.  Barrett,  H.  H.   "Pulse  Compression  Techniques  in  Nuclear  Medicine." 
Procedures  IEEE,  60:723  (1972). 

6.  Rogers,  W.  L. ,  K.  S.  Han,  L.  W.  Jones  and  H.  W.  Beierwaltes. 
"Application  of  a  Fresnel  Zone  Plate  to  Gamma  Ray  Imaging." 
Journal  of  Nuclear  Medicine,  13:612  (1972). 

7.  Sider,  R.  G. ,  E.  S.  Kenney  and  A.  M.  Jacobs.   "Pseudo-Holography 
Using  Scattered  X-Rays."  Nature,  249:546  (1974). 

8.  Rose,  R.  P.   Personal  Communication.   The  Pennsylvania  State 
University,  University  Park,  Pa.  (1974). 

9.  Shulman,  A.  R.   Optical  Data  Processing.   Wiley,  New  York,  p.  387 
(1957). 

10.  Hielstrom,  K.   DROMB,  OS/360  Fortran  IV  Routine  for  Two  Dimensional 
Quadrature.   Argonne  National  Laboratory,  Argonne,  111.   (1966). 

11.  U.S.  Atomic  Energy  Commission.   ANL-5800,  Reactor  Physics  Constants. 
U.S.  Government  Printing  Office,  Washington,  D.C.,  2nd  ed.  (1963). 

12.  Sider,  R.  G.   "The  Use  of  Coded  Aperture  Techniques  in  the  Formation 
of  Images  with  Scattered  Penetrating  Radiation."   M.S.  Thesis, 

The  Pennsylvania  State  University,  University  Park,  Pa.  (1974). 


;r/les/'s 


157339 


Johannesmeyer 

A  contribution  to  th* 
study  of  imaging  with 
scattered  thermal  neu- 
trons. 


3  2768  002  10772  4 
DUDLEY  KNOX  LIBRARY 


